typeset using JPSJ.sty <ver.0.7f> 



Generalization of Gutzwiller Approximation 

Takuya Okabe 

Department of Physics, Kyoto University, Kyoto 606-01 

s ■ 

0> ; (Received January 20, 1997 ) 

pi i' We derive expressions required in generalizing the Gutzwiller approximation to models 
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comprising arbitrarily degenerate localized orbitals. 
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Previously, we generalized the Gutzwiller approximation to degenerate-band models, and 
investigated itinerant ferromagnetism in the 3(i transition metal systems.0) There, we pre- 
sented the resulting formula but omitted its derivation because it was quite complicated.^ 
However, once we know the physical interpretation of the result as explained in ref. |T], it is 



^ . natural to expect that the formula may be derived rather straightforwardly in the way that 
'-^ ■ it reflects the simple interpretation. In fact, as we found this is actually the case, we report 
Q ' it in the following. At the end, we show that the Brinkman-Rice transition occurs only when 



the carrier density equals an integer. 
We investigate the Hamiltonian comprising arbitrary numbers of localized orbitals /; 



H = T + v = Y.t]fi],cu + Y. E c(^n/i,---,g^>fH/i,---,g. (i) 

For simplicity, the hopping part T is assumed to be diagonal with respect to /. In the poten- 
tial part V, indices Im designate orbital and spin of localized states. Operators C'l^\li, ■ ■ ■ ,lp) 
have eigenvalue 1 if only p orbitals h, - ■ ■ ,lp are occupied (the others unoccupied) at site i. 
These operators correspond to di = hiihi^ of the single band Hubbard model. The sum in 
V is taken over a set {li, ■ ■ ■ ,lp}, not to count their combination redundantly. The Pauli 
principle requires Im 7^ for m 7^ n for the orbitals in the braces {/i, ■ ■ ■ , Ip}. C^p\Ii, ■ ■ ■ ,lp) 
is interaction energy of an eigenstate of ■ ■ ■ ,lp). Generalizing the Gutzwiller approx- 
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imation, we can estimate the expectation value of H as 

h^\H\^)=Y^q{l)e^+ Y: C^'\hr--Jp)'y'^'\hr--,lp), (2) 
^ I p 

where ii is defined by 

ei = -{^Q\Y,t'jic]iC,i\^o), (3) 
^ i,j 

and represents the average kinetic energy of the band / for the uncorrelated state I^I/o), for 
which we assume (\l/o|\l/o) = 1- In eqs. (0) and d), L is the total number of lattice sites. 
Parameters u^^^li, ■ ■ ■ ,lp) are the expectation value of z^j-^''(/i, ■ ■ ■ , /p), and take non- negative 
values. In the uncorrelated case, where C^^^ = for any p, these are given by 

z/(^')(Zi,---,g = i.J^)(/i,---,g^ n ^ih) n i^-n{h)), (4) 

i={l,-,p} iy^{l,-,p} 

where 

'^(0 = 7(*olE4cd^o). (5) 

Most difficult and laborious part in applying the Gutzwiller approximation is the calcu- 
lation of the band-width reduction factor q{l) in eq. (|^) as a function of z/(^^(/i, ■ ■ ■ ,Zp). If 
this is achieved, //'•^•'(/i, ■ ■ ■ , /p) are determined so as to minimize {^\H\^). In the variational 
calculation, these parameters z/'-^-' cannot take arbitrary values since they must satisfy two 
relations,0^ i.e., 'conservation of number', 

n(/i) = i/«(/i)+ E E ^^^n^i, ^2, (6) 

p{>2) {/2,-^}(7^/l) 

and 'conservation of probability', 

E E y^'\k.---.ip) = i. (7) 

p {h,-,ip} 

Therefore for fixed n{li) we must regard only i/^^^ with p > 2 as variational parameters. 
A general result for q{l) was obtained for the first time in ref. ^ 

g(^i)= ^nn-.r/ ^ v'^(-)(/i,---,gv'^(^-^)(/2,---,g) , (8) 

and interpreted physically as shown in Fig. 1. Terms in the round bracket squared are inter- 
preted as a sum of products z/(p)\/ of probability amplitudes (a/zTp)"), which 
are involved in the hopping process under consideration. The square is due to contribu- 
tions from two sites i and j. One can make sure that ■ ■ ■ ,lp) defined in eq. meet 
eqs. (I) and (|^), and thus verify that q{l) = 1 for C^^^ = using eq. (0), as expected. To 
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Fig. 1. Examples of the hopping process. A carrier in the orbital li at site i moves to the orbital ly at site 
J- 



derive eq. we must derive relations between u^^^ and the Gutzwiller projection factors 
rjip) 1)^ which are defined in the generahzed Gutzwiller wavefunction, 

1^) = n [i - (i - ^^'^^1' ■ ■ ■ ' ^p)) ^^^'^^1' ■ ■ ■ ' ^ i^o), (9) 

i 

P,{h,---,lp} 

where rj^^^k) = rj^'^^ = 1. 
In the Gutzwiller approximation, the norm of the state (|^) is estimated as 

2pW(«i,-,/p) 



p p(p) 
{h,-,ip} 



iV(L,{iV(/,)},{z7(^)})p(L,{iV(/,)}),(10) 



where 



(/i,---,g^Lz/(p)(/i,---,g, 



N{k)=Ln{h), 



N{L,{N{k)},{u^^^} 



p 

{h,-,lp} 



N{h) 



P{L,{N{k)}) = X{n{k) (l-n(/.)) 



\L-N(h) 



(11) 



These are ereneralization of the results of Vollhardti for the single band model. 

Then, we shall approximate the sum over 1/*^^^ in eq. (10) by a single term which gives the 
largest contribution in the thermodynamic limit L oo; u^^^ {p > 2) is determined by 

d 



from which we obtain a relation between rj^^^li, ■ ■ ■ ,lp) and u^^^li, ■ ■ ■ ,lp), 

2 (z/(o))p-v(^')(/i,---,g 



(12) 



z/«(/i)---/^«(g 



(13) 
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To derive eq. (13) from eqs. (10) and (p!2[), one may note 

dz/«(/,) 



dz/iP)(/i,---,/p) 
dz/(°) 

d.(^)(/„...,g=^-^' 

owing to eqs. (^ and (|^), and 

^logiV! = logiV. (iV^oo) (15) 

As a next step, we note that spatial correlation of various configurations is completely 
neglected in the Gutzwiller approximation. Thus, the factor q can be separated into two 
independent parts, each of which comes from contribution for the creation and annihilation 
operators; 

— f r i r — . (16) 

(^1^) \, % (^ol^o) 

To calculate r for the process shown in Fig. 1, two factors appear besides the Gutzwiller 

h 

parameters r]^^^ ({/i, ■ ■ ■ , Ip}) r/^P"^) ({/g, ■ ■ ■ , Ip}). These are the following ratios, i.e. 
N{L-1, N{h) - 1, Nik) - 1, ■ ■ ■ , Nilp) - 1; ■ ■ ■) _ ^^^'^ ihy'^ (k) ■ ■ ■ ^^'^ 



N{L,N{h),N{k),---,N{lpy,---) (i.(o))P-i 

and 

P(L-l,iV(Zi)-l) 1 



17) 



P(L,iV(/i)) n(/i)' ^ ' 

where dots after the semicolon in the denominator and numerator of the left-hand side of 
eq. (17) represent common parts. We need these ratios since the orbital li at the site i, 
where orbitals h, - ' ' Jp ci-re occupied, should also be occupied before the hopping process. 
Therefore the other N{li) — 1 {i = 1, - ■ ■ ,p) carriers must be on the other L — 1 sites. 
Finally, we obtain the result for rc,^ , 

^«(/i) 



'"'^ n{h) 



1 + 



+ E ./"(0..-.«)^-'>(fc--.«)-'''""^'-''"'<''' 



p>2 



(z/(o))p-i 



(19) 



In the same fashion, in terms of 

iV (L - 1, Njly), N{h) - 1, ■ ■ ■ , Njlp,) - 1; ■ ■ ■) _ u^''^ u'-'Hh) ■ ■ ■ v^^\lp') 
N{L,N{lv),---,N{lp.)----) (z/(o))p'-i 

and 

P(L-l,iV(/i,)) _ 1 
P{L,N{ly)) l-nilvY 



(20) 



(21) 
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one obtains 
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1 + 

, , ,/(i)C7„,^ . . .,/(i)('7 Ai 

■ (22) 



+ E V^''H{h',---,lp'})v^''-'H{l2',---Jr'}) 

p'>2 



(z/(o))p'-i 



Substituting eq. ([T3| ) to eqs. (19) and (22) for /i = ly = /, we finally conclude the result (|D 
for q{l) = r^tr^^. 

Generally in the case I 4^ I', asymmetry arises between r t and r , or, one obtains different 
factors qi'i and qw for (c^ji/Cu) and {cliCji'), respectively. These should be the same because 
of the Hermite character of the Hamiltonian. Thus in this situation, one may instead use 
their average q = {qi>i + qw)/^ for both of {c^ji/Cu) and {c\iCjii). In the same manner, the q 
factors for (c];/c|;) and {cuCji') take different forms. In this case, however, one can use their 
geometrical mean since these factors always appear pairwise as a product in physical terms. 

Notwithstanding general treatment, compared with the method of Vollhardt,!!' our deriva- 
tion presented above is not only simple but transparent with respect to the physical meaning 
of respective terms contributing to the q factor. This is due to observation made at (0); we 
treated the two sites involved separately, while usually these are treated altogether. In par- 
ticular, the square root viAp) appears in eq. (§) in place of r]^P\ eq. (|T3|) , while squared 77^^^ 
there stems from the norm (\E'|\E'), eq. (10). Therefore our derivation validates our previous 
interpretation of q{l) as a sum of products of probability amplitudes V z/(p) . 

In passing, we end this article with a comment on the metal-insulator transition of the 
type noted first by Brinkman and Rice.S* Summing eq. (^ over li, we obtain 

n = EM^i) = E^^'H^i) + E E E ^^'\hj2r--,ir) 
= T.^^'\h)+ E E v^^''\kM.---.ir>) 

h p(>2){h,--,lp} 

= E E P^^''\hMr--,lv)- (23) 
p {h,--,ip} 

If the on-site interaction is strong enough, one may consider only configurations which have 
p— 1, p and p+1 carriers on a single site, where an integer p is selected so that p—1 < n < p+1. 
Then, if we set probabilities of the other configurations as zero, eqs. (23), (^ and (^ become 



n= E {p-i)y^''-'\ii.i2.---.ip-i)+ E P^^^H^i,^2,---,g 

+ E {P + I)i^^'^'\lul2.---Mi)^ (24) 

1= ^^'-'\h,i2,---,ip~i)+ E ^^'\h,i2,---,ip) 
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+ iy^''^'\h,l2,---,M, (25) 

and 



+ E v'z.(P+l)(/i,---,/,+i)v'/^W(/2,---,Wl)) 

Subtracting p times eq. (25) from eq. (24), we get 
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(26) 



n-v=- E v^'-^Kh^h.-'-.K-x)^ E ^(^+'n^i,^2,---,Ui)- (27) 

{'i,-"i'p-i} {^ir"i'p+i} 

For the Brinkman-Rice transition that g(Z), eq. (26), to vanish, the two terms in the right- 
hand side of eq. ( pTl) must vanish at a time. This is because v'^^^ in eq. (26) generally remain 
non-zero to accommodate n particles per site, or to meet eq. (24). (Note that z/^^^ > for 
any p.) Thus we conclude that the Brinkman-Rice transition occurs only when the carrier 
density n equals an integer j9.& 
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